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Abstract 

Pectin is a main component of the plant cell wall and is the most complex family of polysaccharides in nature. Its 
composition is essential for the normal growth and morphology pattern, as demonstrated by pectin-defective mutant 
phenotypes. Besides this basic role in plant physiology, in tomato, pectin structure contributes to very important 
quality traits such as fruit firmness. Sixty-seven different enzymatic activities have been suggested to be required for 
pectin biosynthesis, but only a few genes have been identified and studied so far. This study characterized the tomato 
galacturonosyltransferase (GAUT) family and performed a detailed functional study of the GAUT4 gene. The tomato 
genome harbours all genes orthologous to those described previously in Arabidopsis thaliana, and a transcriptional 
profile revealed that the GAUT4 gene was expressed at higher levels in developing organs. GA(7T4-silenced tomato 
plants exhibited an increment in vegetative biomass associated with palisade parenchyma enlargement. Silenced 
fruits showed an altered pectin composition and accumulated less starch along with a reduced amount of pectin, 
which coincided with an increase in firmness. Moreover, the harvest index was dramatically reduced as a conse- 
quence of the reduction in the fruit weight and number. Altogether, these results suggest that, beyond its role in pec- 
tin biosynthesis, GAUT4 interferes with carbon metabolism, partitioning, and allocation. Hence, this cell-wall-related 
gene seems to be key in determining plant growth and fruit production in tomato. 
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Introduction 

The main role of cell walls is to confer physical strength 
and provide a barrier against the external environment. 
Pectins constitute the ~35% of the primary cell wall in 
eudicot and non-graminaceous monocots, and 2-10% in 
grasses (Mohnen, 2008). Thus, as a main component of the 



cell wall, pectins contribute to these two functions together 
with cellulose and hemicellulose. Pectins include a fam- 
ily of polysaccharides, whose major constituent is galactu- 
ronic acid (GalA) (~70%). Four pectin polysaccharides have 
been described: homogalacturonan (HG), xylogalacturonan 



Abbreviations: AIR, alcohol-insoluble residue; Ara, arabinose; Fuc, fucose; GalA, galacturonic acid; GATL, GAUT-like; GAUT, galacturonosyltransferase; GC-MS, Gas 
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(XGA), rhamnogalacturonan II (RG-II), and rhamnogalac- 
turonan I (RG-I). HG is a linear polymer made of a- 1,4- 
linked GalA and is the most abundant component of pectin 
(~65%). XGA is a modified HG with lateral ramifications of 
xylose attached to the 0-3 position of some GalA residues. 
RG-II is a much more complex structure and makes up ~10% 
of the pectin. In RG-II, the HG backbone is complemented 
by lateral chains composed of 12 different sugars attached by 
22 distinct linkages. Finally, RG-I is the only pectin polysac- 
charide whose backbone is not exclusively made of GalA but 
is composed of a disaccharide (a-l,4-D-GalA-a-l,2-L-Rha) 
repeating unit where the L-rhamnose (Rha) residues may be 
branched with arabinan, galactan, and/or arabinogalactan 
(Harholt et al, 2010). Different models have been proposed 
to explain how these structural elements are combined into a 
macromolecular structure. The model most consistent with 
recent results suggests the existence of an RGI-I backbone 
decorated with side chains composed of the other pectic con- 
stituents (Vincken et al, 2003). 

Galacturonosyltransferases (GAUTs) are among the 
enzymes responsible for pectin biosynthesis and belong to a 
vast group of glycosyltransferases that, together with other 
carbohydrate active enzymes, are catalogued in the CAZy 
database (http://www.cazy.org; Cantarel et al, 2009). More 
specifically, they can be assigned to the glycosyltransferase 
family 8 (GT8) (Harholt et al, 2010; Yin et al, 2010) con- 
sisting of three separate protein classes, class I and class II 
contain mostly eukaryotic proteins, while class III consists 
almost entirely of bacterial proteins. The plant cell-wall- 
related proteins are all located in class I, which includes the 
galacturonosyltransf erase (GAUT) genes and the GAUT-like 
(GATL) genes (Sterling et al, 2006). An extensive phyloge- 
netic analysis revealed that the GT8 family was acquired from 
an ancient cyanobacterium progenitor and further separated 
into subclades, which have undergone functional specializa- 
tion. Moreover, the increase in the numbers of GAUT and 
GATL genes from lower to higher plants suggests that these 
genes have an important role in plant adaptation to living on 
land and adopting an upright growth habit (Yin et al, 2010). 

In Arabidopsis thaliana, 15 GAUT- and 10 GATL-encoding 
genes have been described, and cumulative evidence indicates 
that they are involved in cell-wall pectin and hemicellulose 
biosynthesis (Orfila et al, 2005; Sterling et al, 2006; Caffall 
et al, 2009; Kong et al, 2011). However, redundancy among 
the individual GAUTs and pleiotropic effects complicate the 
assignation of specific functional activities to each protein. 
Given these problems, only four of them have been function- 
ally characterized to date. Bouton et al (2002) reported that 
the A. thaliana mutant quasimodol {qual), which harbours 
a T-DNA insertion in the GAUT8 gene, displays dwarfism 
and reduced cell adhesion. The qua! cell walls contained sig- 
nificantly higher proportions of arabinose (Ara), Rha, and 
fucose (Fuc), and a lower proportion of uronic acid and xylose 
(Leboeuf et al , 2005, Orfila et al. , 2005). The reduction in cell 
adhesion in qual was explained by the low content in calcium- 
dimerized HG, while the alterations of pectin composition 
indicated modifications of other pectic domains like RG-I 
side chains (Leboeuf et al 2005). Moreover, Sterling et al. 



(2006) biochemically characterized the homogalacturonan 
GAUT activity in A. thaliana, and identified GAUT1 protein 
as well as its encoding gene, demonstrating its role in oligoga- 
lacturonide elongation. Recently, it has been demonstrated 
that GAUT1 acts in the Golgi lumen in a complex manner 
together with GAUT7, which anchors the protein complex 
to the Golgi membrane (Atmodjo et al, 2011). Finally, the 
characterization of an irregular xylem8 (ir.\8) dwarf mutant 
allowed the proposition that the wild-type locus product, 
GAUT12, is involved in glucuronoxylan (GX) biosynthesis, a 
hemicellulosic component of the secondary cell wall in dicot 
plants. The mutation led to a 57% reduction in the xylose 
(Xyl) content of GX polysaccharide, and it was proposed 
that GAUT 12 catalyses the addition of an a-D-GalA residue 
to 0-4 of the Xyl residue of the |3-D-Xyl/?-(1^4)-|3-D-Xyk?- 
( 1 — >3)-a-L-Rha/?-( 1 — >2)-a-D-Gal/A-( 1 — >4)-D-Xyl sequence 
present at the reducing end of GXs (Peiia et al, 2007). 

As demonstrated by the characterization of mutant phe- 
notypes, GAUT activity, and consequently pectin composi- 
tion, is essential for normal growth and morphology (Caffall 
et al, 2009). Besides this crucial role, in tomato, the amount 
and composition of pectin and hemicellulose affect fruit firm- 
ness (Giovannoni et al, 1989; Tieman et al, 1992; Brummell 
et al, 1997, I999a,b; Saladie et al, 2007; Chapman et al, 
2012; Lahaye et al, 2012; Lunn et al, 2013). However, no 
GAUT genes have, as yet, been functionally described in this 
species. To gain further insight into the function of pectin 
composition in the whole plant physiology, in this work we 
characterized the GAUT gene family in tomato, focusing on 
a functional study of a putative GAUT-encoding gene asso- 
ciated with the content of free galacturonate in ripe fruits 
(Berimidez et al, 2008). The results demonstrated that this 
gene affects pectin content and solubility, suggesting that dis- 
turbance of cell-wall metabolism causes pleiotropic effects 
leading to significant changes in growth and carbon parti- 
tioning in tomato. 



Materials and methods 

Plant material, growth conditions, and sampling 

Tomato seeds from Solanum lycopersicum L. (cv. M82 and 
Moneymaker) were obtained from the Tomato Genetic Resource 
Center (http://tgrc.ucdavis.edu). Nicotiana benthamiana seeds were 
obtained from Meyer Beck (Berlin). Tomato and tobacco plants 
were grown in 20 and 1 litre pots, respectively. The greenhouse con- 
ditions were 16h/8h photoperiod, 24 ±3 °C, 60% humidity, and 
140 ±40 umol nr 2 s -1 incident irradiance. 

For expression profile analysis, source (the first totally expanded 
leaf, which for our cultivar and growing conditions corresponded to 
the third leaf from the top of the plant) and sink (the first apical leaf 
not fully expanded) leaves were collected from 8-week-old plants (cv. 
M82). Fruit pericarps (without placenta and locule walls) at green, 
mature green, breaker, and ripe stages were harvested 30, 45, 50, and 
60 d after anthesis, respectively. All samples were obtained from six 
independent plants, immediately frozen in liquid N 2 , and stored at 
-80 °C until use. Samples were pooled in three replicates (two inde- 
pendent plants per pool) for further analyses. 

For transgenic plant phenotyping, at least three replicates of the 
selected TO lines were established. After 16 weeks, source leaves 
and ripe fruit pericarps were collected as described above. Fruit 
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biochemical phenotyping was performed exclusively at the ripe stage 
due to the extremely low fruit production of G^(/7¥-silenced plants. 
Six-month-old plants were harvested, and fresh and dry aerial bio- 
mass were determined. These parameters were used to determine 
water content and harvest index according to the following formula: 

Harvest index = (fruit fresh massxlOO)/total aerial biomass. 

For heritability assays, 15 Tl plants for each TO transgenic 
line were evaluated for the presence of the transgene by PCR. As 
expected, the transgene segregated at a 3:1 ratio according to a % 2 
test (P <0.05). GAUT4 expression was reduced by at least 60% for 
all three transgenic lines (Supplementary Fig. S1A at JXB online). 
In agreement with the TO phenotype, Tl silenced plants displayed a 
higher biomass and were taller (Supplementary Fig. SIB; see also 
Fig. 5A). All measurements in the Tl generation were performed in 
triplicate in 8-week-old plants. 



Light microscopy 

Completely expanded leaves were fixed in neutral buffered formalin 
in phosphate buffer (pH 7.0) for 48 h (Lillie, 1965), dehydrated in 
a graded ethanol series, and embedded in plastic resin (Historesin 
Leica; Gerrits, 1991). Cross-sections (10 urn) were cut on a rotating 
microtome (Reichert Jung Autocut 2040; Leica), stained with 0.05% 
(v/v) toluidine blue in 0.1 M sodium acetate (pH 4.7) (O'Brien et at, 
1964), and analysed with an Olympus BX41 light microscope with 
a coupled digital camera. Total, palisade, and sponge parenchyma 
areas were measured using ImageJ software (Rasband, 1997-2012). 



Identification and analyses of GAUT sequences 

The complete sequence of the S. lycopersicum GAUT4 gene, 
Solyc04g0 15270, was obtained from the Solanaceae Genomics 
Network (SGN) database (Bombarely et al, 2011). Untranslated 
regions were determined based on comparison with a reference 
unigene (U320717) using MULTALIN software (http://www- 
archbac.u-psud.fr/genomics/multalin.html; Corpet, 1988). 

Gene family protein sequences were identified in the SGN data- 
base by BLAST search, using the A. thaliana and Oryza sativa ortho- 
logues described by Caffall et al. (2009). GAUT protein sequences 
were aligned using ClustalX (Thompson et al, 1997) and phylo- 
genetic analysis was performed by MEGA5 (Tamura et al, 2011). 
Evolutionary distances were computed using the JTT matrix-based 
method (Jones et al, 1992) and expressed as the number of amino 
acid substitutions per site. The bootstrap values were determined 
from 1000 trials. As outgroups, an A. thaliana (Atlgl9300; Kong 
et al, 2011) and S. lycopersicum (Solyc02g065530) GATL sequence 
were included in the analysis. The identified tomato genes were 
mapped onto the Tomato-EXPEN 2000 genetic map available at the 
SGN database. The genetic positions were obtained by BLASTN 
(Altschul et al, 1990) against the entire Tomato-EXPEN 2000 map 
marker sequence database. Map Chart software 2.2 (Voorrips, 2002) 
was used to construct the graphical representation of the genetic 
map. Protein subcellular localization prediction was performed 
by using TargetP software (Emanuelsson et al, 2007), SherLoc 
(Shatkay et al, 2007), MultiLoc (Hoeglund et al, 2006), Plant- 
mPLoc (Chou and Shen, 2010), and iPSORT (Bannai et al, 2002) 
software. Transmembrane domains were predicted using TMpred 
(http://www.ch.embnet.org/software/TMPRED_form.html). 



RNA isolation and quantitative PCR (qPCR) analyses 

Total RNA from lOOmg of samples was extracted with TRIZOL rea- 
gent (Invitrogen) according to the manufacturer's instructions. DNA 
was removed with 1 ul of amplification-grade DNase (Invitrogen) fol- 
lowing the recommended protocol. cDNA was synthesized from 1 |xg 
of RNA using random primers and a Superscript III kit (Invitrogen), 



and quality was confirmed by PCR using actin intron-fianking spe- 
cific primers (ActinaUp 5'-TGGCATCATACCTTTTACAA-3' and 
ActinaLow 5'-TCCGGGCATCTGAACCTCT-3'). Elongation 
factor 1-a (EFl-a), TIP4, and EXPRESSED were used as house- 
keeping genes. mRNA levels were quantified by real-time qPCR 
using a 7500 Real-Time PCR system (Applied Biosystem), SYBR 
Green Master Mix (Applied Biosystem), and specific primers 
(Supplementary Table SI at JXB online) at final concentrations of 
200nM. PCR conditions were: 95 °C for lOmin, 40 cycles of 95 °C 
for 15 s, primer annealing temperature for lmin, and 72 °C for 30 
s. All reactions were performed with two technical replicates and at 
least three biological replicates. 

Data were analysed with LinRegPCR software (Ruijter et al, 
2009) to obtain C t values and primer efficiency. Relative expression 
and statistics analysis were calculated using fgStatistics software (Di 
Rienzo, 2009). 

Cloning procedures 

Primers were designed using the software Oligo Analyzer 3.1 (http:// 
www.idtdna.com) on the basis of unigene sequences. For the GAUT4 
subcellular localization experiment, a 699 bp fragment spanning the 
first 23 3 aa of the protein was amplified by PCR usingTaq Platinum Pfx 
DNA polymerase (Invitrogen) and the following set of primers FUS- 
GAUT-F ( 5 ' -C ACC ATGA AG ATG AA ACTG AGG A AGCC- 3 ' ) 
and FUS-GAUT-R (5'-TGAACCACCCAATTTGTTG -3'). PCRs 
were conducted in a total volume of 50 ul containing 0.2 mM of each 
dNTP, 0.2 mM of each primer, lx reaction buffer, 1.5mM MgC10 4 , 
50 ng of cDNA, and 2.0 IU of enzyme. The amplification conditions 
were: 94 °C for 3min, 35 cycles of 94 °C for 30 s and 53°C for 30 
s, and 68 °C for 3 min. The amplification product was cloned into 
a pENTR/d-TOPO vector (Invitrogen) and the expression cassette 
was then transferred to the binary vector pK7FWG2 by recombina- 
tion using LR clonase (Invitrogen) resulting in a C-terminal green 
fluorescent protein (GFP) fusion protein (pK7FWG2-GAUTLyc). 

For RNAi silencing, a 161 bp fragment of the GAUT4 gene 
was used to generate a hairpin construct. The fragment was 
amplified with the primers GAUTRNAi-F (5'-CACCCGAGCA 
GCAGTCATACACTAC-3') and GAUTRNAi-R (5'-GTGCAG 
AG A AC ACTTC ATG A ACC AC- 3 ' ) , and Taq DNA polymerase 
(Invitrogen). The amplicon was cloned into the pENTR/d-TOPO 
vector (Invitrogen) and transferred into the pK7GWIWG2(I) 
binary vector using LR clonase, to generate a hairpin construct 
(pK7GWIWG2(I)-GAUTLyc). 

Binary vectors (Karimi et al, 2002) were introduced into 
Agrobacterium tumefaciens strains GV3101 and GV2260 for subcel- 
lular localization and plant stable transformation, respectively. 

N. benthamiana transient transformation and confocal 
microscopy 

The Agrobacterium strains containing pK7FWG2-GAUTLyc, a 
cytosolic control pK7FWG2-HPPDLyc (full-length cDNA frag- 
ment of the hydroxyphenylpyruvate dioxygenase (HPPD) fused to 
GFP; Almeida et al. unpublished results), and a Golgi complex 
control [based on the cytoplasmic tail and transmembrane domain 
of GmManl, soybean a-l,2-mannosidase I fused to yellow fluo- 
rescent protein (YFP); Nelson et al 2007] were grown at 28 °C 
in LB medium with appropriate antibiotics to an optical density 
of 0.6 at 600 nm. The cells were harvested, resuspended in lOmM 
MES with lOOmM acetosyringone (Sigma), incubated for 4h in 
the dark at room temperature, and then infiltrated into leaves of 
6-week-old N. benthamiana plants. After 48 h, the infiltrated tis- 
sues were observed with a confocal laser microscope (Zeiss LSM 
400) under a 63 x water objective. Chlorophyll images were cap- 
tured over 590 nm after excitation at 543 nm, while the GAUT:GFP 
fusion, HPPD:GFP, and Golgi:YFP control were captured over a 
505-550 nm range after excitation at 488 nm with an argon laser 
beam. 
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Plant transformation 

Seedling cotyledons of S. lycopersicum (cv. Moneymaker) were 
used as explants to generate transgenic tomato plants with the 
hairpin construct pK7GWIWG2(I)-GAUTLyc. Agrobacterium- 
mediated transformation was carried out as described previ- 
ously (Nunes-Nesi et al, 2005). The presence of the transgene 
was confirmed by PCR with 35S promoter-specific primers 
(35S-right, 5'-CCCACTATCCTTCGCAAG-3' and 35S-left, 
5 ' -GC AGGTC ACTGG ATTTTGG-3 ' ) . Thirteen independent TO 
transgenic lines were checked for GAUT4 expression level by qPCR. 
Three lines were chosen for further detailed phenotypic charac- 
terization. GAUT4 transcript levels were checked in ripe fruits and 
leaves in at least six biological replicates for each line. 

Photosynthetic parameters 

Photosynthetic parameters, including gas exchange and chlorophyll 
fluorescence emission, were evaluated in 4- 5-week-old plants, main- 
tained at a fixed irradiance of 140 ±40 umol mr 2 s~ l using a Li-Cor 
6400 system (http://www.licor.com/), under different light inten- 
sities (200, 400, 800, and 1200 umol mr 2 s 1 ), 400 ppm C0 2 , and 
a leaf temperature of 25 °C. All measurements were performed in 
dark-adapted plants on the second or third leaflet of the third fully 
expanded leaf. Photosynthetic parameters were estimated according 
to the method of Maxwell and Johnson (2000). Non-photochemical 
quenching (NPQ), photochemical quenching (qP), electron trans- 
port rate (ETR), and reduced plastoquinone accumulation (1 - pQ) 
were estimated using the following equations: 

NPQ = (Fm - Fm')/Fm' 

qP = (Fm'-Fs)/(Fm'-F0) 

ETR = [(Fm' - Fs)/Fm']xflaleaf 

1 - pQ = (Fs - F0')/(Fm' - F0') 

where Fm and Fm' are the maximal chlorophyll fluorescence after 
dark and light adaptation, respectively; Fs is steady-state fluores- 
cence; F0 and F0' are the minimal chlorophyll fluorescence after 
dark and light adaptation, respectively; /is the fraction of absorbed 
quanta that is used by Photosystem II, which is assumed to be 0.5 
for C3 plants; / is the incident photon flux density, and odeaf is leaf 
absorbance, which is assumed to be 0.85 (Sienkiewicz-Porzucek 
et a/., 2010). 

Fruit firmness measurements 

Firmness was measured on similar-sized intact ripe fruits (60 d after 
anthesis). Measurements were performed twice on at least one fruit 
per plant and three plants per line, totalling 44 assessments, accord- 
ing to the method of Smith et al. (2002). Each fruit was tested with 
a cylindrical flat plate (35 mm diameter) with the probe centred over 
different locules at 90° on the equator of the fruit. Force deforma- 
tion curves were recorded using a TA.XT Plus Texture Analyzer 
(Stable Micro Systems) as follows: pre-test speed of 1.5 mm s _1 , test 
speed of 1 mm s~ and a compression distance of 3 mm. Firmness 
was defined as the maximum at a constant deformation (3 mm). 

Soluble sugars and starch quantification 

Freeze-dried samples of leaves and fruits were ground in a ball 
mill (TE 350; Tecnal, Brazil). Ten milligrams of each sample was 
extracted five times with 1.5 ml of 80% (v/v) ethanol at 80 °C for 
20min. The total extraction volume (7.5ml) was dried under vac- 
uum and the samples were resuspended in 1 ml of deionized water. 
To remove pigments, an extraction with 0.5 ml of chloroform was 
performed. A 500 ul aliquot was used for sugar identification and 



quantification by high performance anion exchange chromatogra- 
phy with pulsed amperometric detection (HPAEC/PAD) using a 
Dionex-DX500 system (Dionex, CA, USA) and a CarboPac PA1 
column using isocratic elution of NaOH (200mM). A calibration 
curve was carried out using standard solutions of glucose, fruc- 
tose, sucrose, and raffinose with a concentration range from 50 to 
200 uM (de Souza et al. , 201 3). For starch quantification, the pellets 
obtained after ethanol extraction were washed with water, dried for 
1 h at 60°C and incubated for 30min at 75 °C with 0.5 ml (120U ml" 1 ) 
of thermostable a-amylase (EC 3.2.1.1) from Bacillus licheniformis 
(Megazyme, Bray Co. Wicklow, Ireland) in lOmM MOPS (pH 6.5). 
This procedure was repeated, giving a total of 120 U of a-amylase. 
Samples were cooled to 50 °C and then digested with 15U of 
amyloglucosidase (EC 3.2.1.3) from Aspergillus niger (Megazyme) 
in lOOmM NaAcO buffer (pH 4.5) twice for 30min at 50 °C. The 
enzymatic reaction was stopped by adding 100 ul of 0.8 M percloric 
acid for protein precipitation. After centrifugation for 5min at 10 
OOOg, an aliquot of 50 ul was incubated with 250 ul of Glucose 
PAP Liquiform (Centerlab; Brasil) containing glucose oxidase (~1 1 
000 U mr 1 ), peroxidase (-700 U ml 1 ), 290 mmol l" 1 of D-4-amino- 
antipirine, and 50 mM of phenol (pH 7.5). Using this system, glucose 
oxidase catalyses glucose oxidation producing H 2 0 2 that reacts with 
d -4-aminoantipirine and phenol under the action of peroxidase. 
These reactions produce a red antipirilquinonimine whose colour 
intensity is proportional to the glucose concentration. After 15min 
at 30 °C, absorbance was measured on a spectrophotometer coupled 
to an ELISA reader at 490 nm. 



Cell -wall composition 

Freeze-dried samples of leaves and fruits (250 mg) were milled 
for further extraction. After extraction with 80% (v/v) ethanol at 
80 °C for five times to remove soluble sugars, 20 ml of 90% (v/v) 
DMSO was added to the pellets and the suspensions were stirred for 
24 h to remove starch, followed by two additions of 20 ml of 90% 
(v/v) DMSO stirred for 3h. The supernatants were discarded and 
the pellets were washed six times with distilled water, freeze dried, 
and weighed. The remained cell-wall material was termed alcohol- 
insoluble residue (AIR). AIR was extracted three times in 20ml of 
aqueous 0.5% (m/v) ammonium oxalate (pH 7) at 80 °C for 1 h each 
with stirring. The supernatants (ammonium oxalate-soluble frac- 
tion) were combined and dialysed (membrane cut off 12 400) for 
1 2 h against running tap water and then dialysed for 8 h with distilled 
water, freeze dried, and weighed. To the remaining cell-wall material 
from leaves, 20 ml of 0.1 M NaOH supplemented with 3mg mF 1 of 
NaBH 4 was added and the suspension was stirred for 1 h. This pro- 
cedure was repeated twice. The supernatants (NaOH-soluble frac- 
tion) were neutralized with glacial acetic acid and dialysed for 12 h 
with running tap water followed by distilled water for lOh. The sam- 
ples were freeze dried and weighed. For both fractions, 2mg aliquots 
were hydrolysed by adding 100 ul of 72% H 2 S0 4 . After 30min at 
30 °C, the solution was diluted to 3% H 2 S0 4 and autoclaved for 1 h 
at 121 °C (Saeman et al, 1945). After deionization through cation- 
and anion-exchange columns (Dowex), the neutral monosaccharide 
composition was determined by HPAEC/ PAD in a Dionex DX-500 
system using a CarboPac PA1 column. Monosaccharides were 
eluted in water using a post-column addition of 500 mM NaOH for 
sugar detection over 40min and a flow rate of 1 ml mhr 1 . 

Uronic acid determination 

The amount of uronic acid was quantified in the ammonium oxa- 
late-soluble fraction essentially as described by Filisetti-Cozzi and 
Carpi ta (1991). Five milligrams of the freeze-dried material was 
hydrolysed with 1 ml of H 2 S0 4 . After 5min in an ice bath, the pro- 
cedure was repeated, giving a total of 2 ml of acid. Subsequently, 
0.5 ml of water was added and after 5min on ice, the procedure 
was repeated. Finally, the mixture was diluted to 10 ml with dis- 
tilled water. In order to perform the colorimetric assay, 40 ul of 4 M 
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K-sulfamic acid/potassium sulfamate (pH 1.6) and 2.4ml of H 2 S0 4 
containing borate were added to 400 of the diluted solution. The 
mixture was boiled for 20min. After the addition of 80 ill of 0.15% 
(m/v) m-hydroxybiphenyl in NaOH 0.5% (m/v), the absorbance 
was read at 525 nm. For each sample, measurements were done in 
duplicate. 



Ascorbic acid 

One hundred milligrams of fresh leaf or fruit tissue was extracted 
with 1 ml of 3% (v/v) trifluoracetic acid (TFA). Five hundred micro- 
litres of the supernatant was applied to an Extract-clean CI 8 col- 
umn (500mg in 8.0ml, 50 urn average particle size, and 60 A pore 
size; Alltech Associates, USA) equilibrated with lOOmM phosphate 
buffer (pH 7), and column was then eluted with 1.5ml of lOOmM 
phosphate buffer (pH 7). To 500 ill of the eluate, 500 ill of lOOmM 
K 2 HP0 4 buffer (pH 8.5) was added. To quantify the total ascor- 
bic acid, 15 ul of 0.1 M dithiothreitol and 125 ul of 3% (v/v) TFA 
were added to 500 ul of extraction mix. To measure the reduced 
ascorbic acid, 15 ul of lOOmM phosphate buffer (pH 7) and 125 ul 
of 3% (v/v) TFA were added to 500 ul of the extraction mix. Both 
samples were incubated for 3min at room temperature. Ascorbic 
acid quantification was performed by HPAEC/PAD using a silica- 
based reversed-phase C18 column (particle size 5 pm, 150 x 4.6mm, 
HL90-5s, Bio-Sil; Bio-Rad, Munich, Germany). The mobile phase 
consisted of a KH 2 P0 4 buffer (lOOmM) at pH 3.0 (with phosphoric 
acid) and was delivered isocratically at a flow rate of 0.5 ml mur 1 . 
Ascorbic acid resulted in a peak at 3.5min. Quantification was done 
based on standard solutions of reduced and oxidized ascorbic acid. 



Gas chromatography '/mass spectrometry (GC-MS) 

For metabolic profiling, frozen leaves (~100mg) were extracted 
in 1.4 ml of methanol, as described by Roessner et al. (2001) with 
the modifications proposed by Lisec et al (2006); 60 ul of internal 
standard (0.2mg ml of ribitol) was added for quantification. In 
the case of ripe fruits (~250mg), the reagent volumes were adjusted 
using 2 ml of methanol and 120 ul of ribitol, as described by 
Osorio et al (2012). In both cases, the mixture was extracted for 
15min at 70 °C, mixed vigorously with 1 volume of water, centri- 
fuged at 2200g-, and subsequently vacuum dried. The residue was 
redissolved and derivatized for 120min at 37 °C (in 60 ill of 30mg 
ml -1 methoxyamine hydrochloride in pyridine), followed by a 30min 
treatment at 37 °C with 120 ill of 7V-methyl-iV-[trimethylsilyl]trif- 
luoroacetamide. Sample volumes of 1 ul were then injected in split- 
less and split modes, using a hot-needle technique. The GC-MS 
system was composed of an AS 2000 autosampler, a GC 6890N 
gas chromatographer (Agilent Technologies, USA), and a Pegasus 
III time-of-flight mass spectrometer (LECO Instruments, USA). 
Chromatography was performed with an MDN-35 capillary col- 
umn of 30 m length, with a 0.32 mm inner diameter and 0.25 urn 
film thickness (Macherey-Nagel). Samples were injected at 230 °C in 
splitless mode with helium carrier gas flow set at 2 ml min -1 . The flow 
rate was kept constant with electronic pressure control enabled. The 
temperature program was isothermal for 2min at 85 °C, followed by 
a 15 °C min -1 ramp to 300 °C. The transfer line temperature was set 
at 250 °C. For detection, the temperature of the ion source was set at 
250 °C (match transfer line condition). Mass spectra were recorded 
at 20 scans s with a scanning range of 70-660 mlz. The remaining 
monitored chromatography time proceeded with the filament turned 
off. The filament bias current was -70 V and the detector voltage 
was 1525 V. Both chromatograms and mass spectra were evaluated 
using ChromaTOF chromatography processing and mass spectral 
deconvolution software, version 3.00 (LECO Instruments, USA). 
Identification and quantification were performed with TagFinder 
4.0 software and the mass spectra were cross-referenced with those 
in the Golm Metabolome Database (Kopka et al, 2005; Schauer 
et al, 2005). Three to six biological replicates were used for this 
analysis. 



Data analyses 

Differences in phenotypic parameters were analysed by Infostat 
software (Di Rienzo et al., 2011). When the data set showed homo- 
scedasticity, an analysis of variance followed by a Tukey or Dunnett 
test (P <0.05) was used to compare transgenic lines against wild- 
type controls. Due to lack of homoscedasticity in some cases, a 
non-parametric comparison was also performed by applying the 
Kruskal-Wallis test (P <0.05). 

Results 

GAUT gene family diversity in S. lycopersicum 

On the basis of a metabolic quantitative trait loci (QTL) 
analysis performed in a Solatium pennellii introgression line 
collection (Schauer et al., 2006), we previously identified a 
putative GAUT-encoding gene co-localizing with a QTL for 
fruit galacturonate content (Bermudez et al, 2008). As no 
GAUT genes have been described in tomato and they belong 
to a large gene family, we first performed a complete phylo- 
genetic analysis. The recently published tomato genome (The 
Tomato Genome Consortium, 2012) was screened to identify 
all GAUT sequences using A. thaliana and O. sativa ortholo- 
gous genes as baits (Caffall et al, 2009). This search retrieved 
17 tomato genes. In order to predict the subcellular localiza- 
tion of the encoded tomato proteins, five different programs 
were tested using At3g61130 (AtGAUTl) and At2g38650 
(AtGAUT7) as controls, which have been demonstrated 
experimentally to localize in the Golgi apparatus (Atmodjo 
et al, 2011). Interestingly, Sherloc, which predicts using 
sequence and text-based features (Shatkay et al, 2007), and 
MultiLoc, which is based on N-terminal targeting sequences 
(Hoeglund et al, 2006), were those that accurately predicted 
both proteins to be located within the Golgi. With the excep- 
tion of S1GAUT2, for which both programs indicated mito- 
chondrial localization, and S1GAUT7, which according to 
MultiLoc is targeted to the endoplasmic reticulum, all tomato 
GAUTs seemed to be directed to the Golgi (Table 1). 

For each A. thaliana GAUT gene, a corresponding tomato 
orthologue was identified, with the exception of the GAUT12 
and GAUT 13 clades, which were rather represented by two 
paralogues in S. lycopersicum. The tree topology demon- 
strated that the GAUT identified by Bermudez et al. (2008) is 
a GAUT4. The clades for GAUT3, -10, -11, and -15 revealed 
a single copy in the eudicot and monocot species analysed, 
and the topology was according to the phylogenetic relation- 
ships between them. The GAUT1/2 and GAUT8/9 clades 
seem to have diverged following the mono/eudicot split, 
undergoing duplication in the monocot lineage as evidenced 
by the presence of three rice paralogues. The GAUT4 and 
GAUT7 clades displayed a single copy in eudicots while being 
represented by three and five paralogues in rice, respectively. 
Finally, clades GAUT5/6 and GAUT13/14 exhibited a similar 
pattern wherein, after the divergence of tomato and A. thali- 
ana, gene duplication in the eudicots was followed by a subse- 
quent gene divergence (Fig. 1A). The 17 tomato GAUT genes 
were mapped onto the Tomato-EXPEN 2000 map housed at 
the SGN website (http://solgenomics.net/) and were found to 
be evenly distributed among chromosomes 1, 2, 3, 4, 6, 7, 10, 



2454 | De Godoy et al. 

Table 1. Identification of tomato genes encoding for galacturonosyltranferases. 



enzyme 


VI f /i n //in i Ia^i io 
M. mallana lOCUb 


lOITldlO 


Subcellular localization 




lumdlO 




o n ro mobome 




(no. amino acids) 


unigene^ 


TargetP c 


SherLoc d 


MultiLoc 


locus 


marker 0 


position' 7 (cM) 


GAUT1 


At3g61130(673) 


SGN-U565384 


nd 


G 


G 


Solyc01g093970 


P51 


1 (69 cM) 


GAUT2 


At2g46480 (528) 


SGN-U598345 


nd 


M 


M 


Solyc10g017600 


T1720 


10(32.5 cM) 


GAUT3 


At4g38270 (680) 


SGN-U567225 


SP 


G 


G 


Solyc01g1 12210 


C2_At4g38240 


1 (165 cM) 


GAUT4 


At5g47780(616) 


SGN-U575018 


SP 


G 


G 


Solyc04g015270 


T0891 


4 (53 cM) 


GAUT5 


At2g30575 (610) 


SGN-U574664 


SP 


G 


G 


Solyc07g005360 


T1112 


7 (0.4 CM) 


GAUT6 


At1g06780(589) 


SGN-U569431 


SP 


G 


G 


Solyc12g010200 


CLET-8-K4 


12 (41 cM) 


GAUT7 


At2g38650(619) 


SGN-U574191 


nd 


G 


ER 


Solyc10g074650 


T0283 


1 0 (36 CM) 


GAUT8 


At3g25140(559) 


SGN-U565194 


nd 


G 


G 


Solyc06g083310 


U146140 


6 (97.2 cM) 


GAUT9 


At3g02350(561) 


SGN-U563256 


M 


G 


G 


Solyc02g089440 


CLES-3-G11 


2(1 1 1 CM) 


GAUT10 


At2g2081 0(466) 


SGN-U574750 


M 


G 


G 


Solyc04g064490 


T0877 


4 (67 cM) 


GAUT1 1 


At1g1 8580 (332) 


SGN-U565457 


M 


G 


G 


Solyc03g1 14810 


CLED-19-N16 


3(107cM) 


GAUT12 


At5g54690 (535) 


SGN-U601760 


nd 


G 


G 


Solyc07g064050 


TG438 


7 (73 cM) 






SGN-U599179 


nd 


G 


G 


Solyc10g006500 


T0418 


10(1 cM) 


GAUT13 


At3g01040(532) 


SGN-U582378 


nd 


G 


G 


Solyc02g067060 


T0869 


2 (21 CM) 






SGN-U563601 


nd 


G 


G 


Solyc03g01 3630 


C2_At3g01060 


3 (72.7 CM) 


GAUT1 4 


At5g 15470 (432) 


SGN-U573429 


nd 


G 


G 


Solyc02g088630 


T1 480 


2(106 CM) 


GAUT15 


At3g58790 (540) 


SGN-U573007 


nd 


G 


G 


Solyc07g055930 


C2_At3g58790 


7 (44.6 CM) 



a Enzyme name abbreviation. 

b Tomato unigene number according to SGN (http://solgenomics.net/). 

c Subcellular localization prediction according to TargetP 1 .1 software (http://www.cbs.dtu.dk/services/TargetP/; Emanuelsson etal., 2007). 
SP, secretory pathway; M, mitochondrion; nd: not determined because output score was below the requested cut-off for localization prediction. 
d Subcellular localization prediction according to SherLoc (Shatkay etal., 2007). G, Golgi; M, mitochondria. 

e Subcellular localization prediction according to MultiLoc (Hoeglund et al., 2006). G, Golgi; M, mitochondria; ER, endoplasmic reticulum. 
' S. lycopersicum locus according to SGN (http://solgenomics.net/). 
9 Closest mapped markers. 

h Chromosome and genetic position in Tomato-EXPEN 2000 v52. 



and 12, with no obvious association between genomic locali- 
zation and phylogenetic distribution (Fig. IB). 

In order to investigate organ expression specificity, a qPCR 
analysis was performed in leaves and ripe fruits. GAUT8, -12, 
-13, -14, and -15 showed higher expression in vegetative tis- 
sue in comparison with ripe fruits, whereas the rest of the 
paralogues did not display differential expression between 
the tested tissues (Fig. 2A). Furthermore, a broader devel- 
opmental expression profile was carried out for the GAUT4 
gene, which showed it to be ubiquitously expressed. However, 
increased levels of mRNA were detected in growing tissues 
such as sink leaves, mature green, and green fruits, although 
in the latter the increase was not statistically significant 
(Fig. 2B). 

GAUT4 structural characterization 

The identification of the genomic clone from S. lycopersicum 
and its comparison with the corresponding mRNA sequence 
allowed determination of the GAUT4 gene structure. GAUT4 
was found to be composed of 12 exons with three being non- 
coding exons. The coding region spanned 2040 nt resulting in 
a protein of 679 aa (Fig. 3A). The protein topology displayed 
both transmembrane and GT8 domains. Within the latter, 
the characteristic GT8 motif DxD was present. Moreover, 
the GAUT family-specific motif described previously by 
Sterling et al. (2006) was also identified, albeit showing an 
amino acid substitution with respect to the Arabidopsis/hce 



consensus sequences. In silico analysis of the protein primary 
structure suggested that GAUT4 is targeted to the Golgi 
(Table 1) and probably remains anchored to the membrane, 
as evidenced by the presence of a transmembrane domain 
(Fig. 3B). This result is in agreement with the localization of 
the GAUT4:GFP fusion protein for which the punctate fluo- 
rescence signal in the cytoplasm of mesophyll cells matched 
that which would be anticipated for proteins targeted to this 
subcellular compartment (Fig. 3C). 

Generation and primary characterization of 
GAUT4-silenced plants 

With the aim of gaining further knowledge about GAUT4- 
specific function, S. lycopersicum transgenic lines expressing a 
GAUT4 hairpin construct were generated. The level of silencing 
was assessed by qPCR in source leaves and ripe fruits, and three 
lines showing significantly reduced levels of GAUT4 mRNA 
(RNAi2, RNAi40, and RNAi42) were propagated for further 
phenotypic characterization (Fig. 4). As functional redun- 
dancy has been hypothesized within the GT family (Mohnen, 
2008), the expression levels of the other four uncharacterized 
GAUT genes (GAUT3, -7, -10, and -7J) were investigated in 
the leaves and fruits of the transgenic lines. These genes were 
chosen because they belong to different clades of the GAUT 
phylogeny spanning all the gene family diversity. Neither co- 
silencing effects nor compensatory upregulation were detected 
for any of the tested genes (data not shown). 
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Fig. 1 . The GAUT protein family and genomic localization of GAUT genes in tomato. (A) Phylogenetic analysis of the GAUT family in 
A. thaliana (At), O. sativa (Os), and S. lycopersicum (SI). A. thaliana and O. sativa sequences were obtained from Caffall etal. (2009). The 
tomato orthologues were identified in this work. The tree was rooted with AtGATL and SIGATL, which are GATL proteins from A. thaliana and 
S. lycopersicum, respectively. Bootstrap values greater than 50% are indicated on the branches. The black dot identifies the gene described 
by Bermudez ef al. (2008). (B) The 1 7 identified genes were localized in the Tomato-EXPEN 2000 genetic map available at the Solanaceae 
Genomics Network (http://solgenomics.net/index.pl). Markers and GAUT genes are indicated on the left side of the chromosomes. GAUT 
genes are highlighted in the same colour pattern as their branches on the phylogenetic tree. Genetic distances are indicated on the left border. 
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(A) 




Fig. 2. Expression profiles of GAUT genes in tomato. (A) Gene expression profiles of 1 4 members of the GAUT gene family were 
measured by qPCR in samples from source leaves and ripe fruits (n=3). Relative gene expression levels (log 2 ) are shown using a colour 
scale where green and red indicate decreased and increased transcript levels in ripe fruits with respect to source leaves (black boxes), 
respectively. Grey boxes indicate that the gene was not detected (nd). The heat map was constructed using MeV software (Saeed 
etal., 2003). GAUT12 and GAUT13 represent the expression of both paralogue pairs: GAUT12-1 and GAUT12-2, and GAUT13-1, 
and GAUT13-2. (B) GAUT4 gene expression profile in source and sink leaves, and in green, mature green, breaker, and ripe fruits. The 
data indicate relative expression normalized to source leaf value (black bar). All data are means from at least three biological replicates. 
Significant statistic differences from source leaves are indicated: *P <0.05. 



Carbon assimilation and other photosynthetic parameters 
were evaluated at the vegetative stage of the plants and no sig- 
nificant differences were observed between silenced and control 
plants (Supplementary Table S2 at JXB online). Despite these 
results, G^t/T^-silenced plants showed higher levels of veg- 
etative growth, as revealed by the size of 8-week-old Tl plants 
(Fig. 5A and Supplementary Fig. SI). The biomass increment 
(Fig. 5B) was mainly due to an increase in the leaf water con- 
tent, which increased ~4% in the silenced plants (Fig. 5D). 
To gain insights about the origin of the biomass increment, 
morphometric analyses were performed from source leaf 



cross-sections of the Tl plants. The palisade parenchyma was 
larger in transgenic plants compared with wild-type controls 
(Supplementary Fig. S2A, B). Furthermore, the cell density 
was also higher in this tissue (Supplementary Fig. S2C) in 
accordance with the differences in growth parameters. 

On assessing fruit characteristics, while flowering time 
was unaltered, the number (Fig. 5E) and weight (Fig. 5F) of 
fruits, as well as the harvest index (Fig. 5G) decreased sig- 
nificantly in the transgenic plants. However, transgenic ripe 
fruits were more resistant to pressure, as demonstrated by the 
flat plate compression test (Fig. 5H). Due to insufficient fruit 
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Fig. 3. Structure of the tomato GAUT4 gene and protein localization. (A) Structure of the tomato GAUT4 gene. Exons are indicated 
as E. Coding and non-coding exons are in black and grey, respectively. Lines represent introns. (B) Schematic representation of the 
GAUT4 protein structure exhibiting the GT8 domain (Glyco_transf 8) as a grey box with the family characteristic motifs, GT8 and GAUT. 
The dark line indicates the position of the transmembrane domain. (C) Confocal imaging of N. benthamiana mesophyll cells expressing 
GAUT4GFP, HPPDGFP (a cytosolic control), and Golgi:YFP (a Golgi complex control) fusion proteins. GFP, chlorophyll fluorescence, 
and merged signals are indicated above the panels. Bars, 25 u.m for GAUT4GFP and HPPDGFP fusions; 50 |im for Golgi:YFP control. 




Fig. 4. GAUT4 expression levels in silenced transgenic plants. GAUT4 expression was measured by qPCR in samples from source 
leaves (A) and ripe fruits (B). Data indicate relative expression means from at least three biological replicates normalized to the wild-type 
(WT) value. Significant statistic differences from the wild type are indicated: *P <0.05. 
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Fig. 5. Phenotypic characterization of the G/AL/7"4-silenced plants. (A) Height of wild-type (WT) and T1 silenced plants at 8 weeks old 
from the base to the apex. (B-D) Six-month-old TO plants were harvested and fresh (B) and dry (C) aerial biomass, and leaf water content 
(D) were calculated. (E-G) For fruits, the total number of fruits (E) and total fruit weight (F) at harvest time were measured, and the harvest 
index (G) was calculated. (H) Firmness was measured for intact fruits at the ripe stage. Values are means from at least three biological 
replicates for (A)-(G), and from at least five fruits per line for (H). Significant statistic differences from the wild type are indicated: *P <0.05. 



production of G^4C/7V-silenced plants, further biochemical 
characterization was performed exclusively at the ripe stage, 
impeding evaluation across fruit development. 

Biochemical characterization of GAUT4-silenced plants 

As GAUT4 was expected to be involved in pectin synthesis, a 
cell-wall fractionation was performed in order to evaluate the 
impact of this gene silencing on the cell-wall composition. 



From leaf tissue, two cell-wall fractions were evaluated in 
which all pectic polysaccharides were supposedly extracted: 
the chelator-soluble fraction, extracted with ammonium oxa- 
late; and the alkali-soluble fraction, extracted with 0.1 M 
NaOH. In the latter, hemicellulose traces could also be pre- 
sent. After the ammonium oxalate extraction, an unwork- 
able amount of material remained from fruit pericarp, and 
thus only the chelator fraction was assessed for this tissue. In 
leaves of silenced plants, the cell wall was basically unaltered, 
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both in amount (Fig. 6) and composition (Table 2). By con- 
trast, transgenic fruits contained less ammonium oxalate- 
soluble pectin (Fig. 6) with reduced content of Rha, Ara, 
and galactose (Gal) (Table 2). For interpretation of the 
occurrence and properties of cell-wall polysaccharides from 
sugar composition, three sugar ratios were analysed: a meas- 
ure for the linearity of pectin, the contribution of RG-I to 
the entire pectin population, and the extent of branching of 
RG-I (Houben et al. , 201 1). The ratios suggested that GAUT4 
silencing affects RG-I accumulation, specially based on Rha/ 
uronic acids ratio, with a clear stronger impact on fruits. 
Interestingly, in both analysed organs, RG-I branching in 
transgenic plants was increased by 26% (Table 3). 

To identify whether the reduction in GAUT4 expression 
could be associated with a shift of plant metabolism in gen- 
eral, a detailed evaluation of primary metabolism was car- 
ried out by GC-MS and HPLC. The resultant metabolic 
profile showed that GAUT4 silencing did not produce mas- 
sive biochemical changes but rather only specific alterations 
in a handful of metabolites (Supplementary Table S3 at 
JXB online). It is worth noting that free GalA was slightly 
increased in leaves, while the increment was more significant 
in fruits in agreement with the lower pectin amount observed 
in the latter organ (Fig. 7A, B). Raffinose was the only altered 
soluble sugar displaying consistently decreased levels in the 
leaves of the transgenic plants (Fig. 7C). 

Finally, in an aim to understand better the relationship 
between the increase in vegetative biomass and the strong 
reduction in harvest index, the starch contents were enzy- 
matically quantified. Whereas in leaves no significant changes 
were observed, in fruits the accumulation of this polymer was 



reduced, suggesting a shift in carbon partitioning between 
source and sink tissues (Fig. 8). 

Discussion 

GAUT genes encode a set of enzymes belonging to the CAZy 
GT8 family of glycosyltransferases that are associated with 
pectin and hemicellulose biosynthesis (Cantarel et al, 2009). 
In the present study, we performed genomic identification of 
the GAUT gene family in S. lycopersicum and functional char- 
acterization of the GAUT4 member. After the identification 
of all GAUT genes encoded in the tomato genome, the phy- 
logenetic analysis revealed that the tomato genome harbours 
all the genes orthologous to those previously described for 
A. thaliana. Caffall et al. (2009) reported the existence of two 
copies of GAUT 12 in poplar, while no orthologues were iden- 
tified in rice. In tomato, our survey retrieved also two GAUT12 
genes, re-enforcing the hypothesis that GAUT 12 has a special- 
ized role in the synthesis of secondary wall GX of dicot walls, 
as proposed by Persson et al. (2007). The tree topology also 
suggested that GAUT 12 constitutes a monophyletic group 
with GAUT13 and GAUT14. This finding is in discrepancy 
with the phylogeny proposed by Caffall et al. (2009), in which 
GAUT15 appeared more closely related to the GAUT13/14 
clade. Three pairs of recently duplicated genes that originated 
after the split between A. thaliana and S. lycopersicum were 
identified: SIGAUT5 and SIGAUT6, and SIGAUT12-1 and 
SIGAUT12-2, as well as SIGAUT13-1 and SIGAUT13-2 and 
S1GAUT14. The 17 tomato GAUTs were mapped, and inter- 
estingly those pairs of genes that were closely linked according 
to the phylogenetic analysis were distantly located within the 
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Fig. 6. Cell-wall polysaccharide content in GAL/7"4-silenced plants. (A, B) The alcohol-insoluble fraction (AIR) containing the cell wall was 
extracted from leaves (A) and fruits (B) and expressed relative to dry material. (C, D) AIR was fractioned by subsequent extractions with 
ammonium oxalate and NaOH from source leaves (C) and only with ammonium oxalate from ripe fruits (D). The mass for each fraction 
was expressed relative to AIR. Significant statistic differences from the wild type for any given fraction are indicated: *P<0.05. 
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genome, suggesting that they have arisen from a segmental 
duplication followed by a genomic rearrangement of the 
tomato genome (Song et al., 2012). A similar observation 
has also been reported for tomato tocopherol biosynthetic 
genes (Almeida et al, 2011). Moreover, Kong et al. (2011) 
determined the relationship between GATL paralogous genes 
also based on the segmental duplication history of the cor- 
responding genomic regions in A. thaliana (Blanc and Wolfe, 
2004). 

Expression profiles suggested that, although cell-wall bio- 
synthesis/degradation is a ubiquitous metabolic process in 
plant cells, there is a certain degree of organ/tissue specificity 
in GAUT activity probably associated with differential pectin 
composition. Tissue specificity has previously been demon- 
strated for GAUT genes. In this sense, by using in situ PCR, 
Orfila et al. (2005) demonstrated that GAUT8 from A. thali- 
ana is expressed in discrete regions around the vascular bun- 
dles. Similarly, GAUT 12 was also expressed in vascular tissues 
in A. thaliana, as revealed by promotengiw fusion studies 
(Persson et al, 2007). In tomato, GAUT8 and GAUT12 are 
more highly expressed in leaves than in fruits. As vascular tis- 
sue density is much more abundant in leaves than in the fruit 
pericarps, our results suggest a similar expression pattern to 
that shown for the A. thaliana orthologous loci. Interestingly, 
GAUT13, GAUT14, and GAUT15, which cluster together 
with GAUT12, had a similar expression pattern, possibly indi- 
cating similar function as proposed by Caffall et al. (2009). In 
particular, GAUT4 was expressed ubiquitously, as was also 
observed in A. thaliana. However, higher expression levels 
were documented in growing stages, pinpointing its crucial 
role in normal development, as indicated by the lack of abil- 
ity to rescue the corresponding A. thaliana mutant (Caffall 
et al, 2009). 

Pectic polysaccharides are synthesized in the Golgi appa- 
ratus of plant cells, sorted to vesicular compartments, and 
subsequently secreted to the apoplastic space (Caffall and 
Mohnen, 2009). Phylogenetic analysis, structural charac- 
terization, and subcellular localization experiments dem- 
onstrated that the gene under study is indeed the tomato 
GAUT4 gene. Moreover, the results are in close agreement 
with the assumed model for biosynthesis of pectin/hemicel- 
lulose cell-wall polymers in which enzyme active sites face the 
Golgi lumen, being anchored, directly or indirectly, by trans- 
membrane domains (Pena et al, 2007; Caffall and Mohnen, 
2009; Atmodjo et al, 2011). 

Pectin comprises as much as 30% of dicot cell walls (Ridley 
et al, 2001). Two-thirds of the worldwide tomato production 
concerns processing tomato, with modifications in the structure 
and composition of pectin being directly related to alterations 
in juice texture and extractability (Sila et al, 2009). Moreover, 
intact fresh fruit texture is a quality index for consumption in 
natura (Barrett et al. , 1998; Lahaye et al. , 2012). To gain further 
insights about pectin metabolism and particularly the role of 
the GAUT4 gene, tomato transgenic plants expressing a hairpin 
RNA construct were obtained. The ripe fruits of silenced lines 
displayed significantly lower levels of ammonium oxalate-sol- 
uble pectin and the pectin contained reduced contents of Rha, 
Gal, and Ara. These data, together with the reduction of 70% 
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Table 3. Sugar ratios for the different polysaccharide fractions of transgenic lines. Ratios were calculated based on Table 2 data. 



Linearity of pectin Contribution of RG-I to Branching of RG-I 

pectin population 

Uronic acids/ Rha/uronic acids (Ara+Gal)/Rha 
(Xyl+Gal+Rha+Ara+Fuc) 

Leaves RNAi plants 91.2 0.0026 3.1 

WT plants 91.2 0.0031 2.3 

Fruits RNAi plants 74.3 0.0006 18.6 

WT plants 28.5 0.0020 13.7 
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Fig. 7. Soluble metabolite content in G/4L/T4-silenced plants. (A, B) GalA was quantified by GC-MS from source leaves (A) and ripe 
fruits (B) and data were normalized to wild type (WT) and expressed as relative content. (C) Free raffinose content in source leaves was 
measured by HPLC and expressed relative to dry weight. This monosaccharide was not detected in ripe fruits in agreement with the 
decrease across fruit development reported by Oms-Oliu et al. (201 1). All values indicate means from at least three biological replicates 
and significant statistic differences from the wild type are indicated: *P<0.05. 
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Fig. 8. Starch content in G4L/7"4-silenced plants. Starch was quantified enzymatically in source leaves (A) and ripe fruits (B). Values are 
expressed relative to dry weight. Data indicate means from at least three biological replicates. Significant statistic differences from the 
wild type are indicated: *P<0.05. 



2462 | De Godoy et al. 



in Rha/uronic acids ratio, suggest that GAUT4 plays a role in 
RG-I backbone synthesis. The absence of the observed effect in 
leaves could be explained by the compensation of other para- 
logues that are also expressed in leaves but not in fruits ( i.e. 
GAUT2, -5, -13, -14, or -75). These other enzymes could either 
have a similar biochemical function and/or culd balancing the 
lack of RG-I producing other polymers or altering branch- 
ing pattern, leading to slight alterations in sugar ratios. In this 
sense, in both analysed organs, an increased level (26%) of 
RG-I branching was observed. Transgenic potato plants with 
a reduced proportion of RG-I displayed an increase in uronic 
acid content also exposing pectin biosynthesis compensation 
(Oomene/ al, 2002). 

In contrast to previous reports on GAUT-deficient plants, 
which displayed dwarfism (Bouton et al. 2002; Pena et al. 
2007; Persson et al. 2007), GAUT4 silencing resulted in 
taller plants with higher vegetative biomass. The morpho- 
anatomical analysis of the leaves highlighted the increased 
water content and fresh mass observed in transgenic plants. 
The reduction in GAUT4 activity might have resulted in a 
looser cell wall; thus, in order to keep the pressure potential 
and consequently the water potential constant, the compact 
palisade parenchyma became enlarged and increased the cell 
density. This is in line with evidence indicating that hydraulics 
play an essential but passive role in controlling cell growth 
(Szymanski and Cosgrove, 2009; Alim et al, 2012). 

Intriguingly, the silencing of GAUT4 resulted in a decreased 
provision of carbon to sink tissues. This was evidenced by 
the decreased harvest index, mostly as a consequence of the 
reduction in total fruit weight and number associated with 
a lower fruit starch accumulation. Recently, Wormit et al. 
(2012) demonstrated that carbohydrate metabolism is respon- 
sive to changes in cellulose biosynthesis via an osmosensitive 
mechanism. All phenotypic data presented here including 
growth parameters, biochemical profile, and histological anal- 
yses suggested that the pleiotropic effects caused by GAUT4 
silencing are the consequence of the same metabolic cross- 
talk, and they allow extension of the proposed mechanism to 
cell-wall biosynthesis. 

Regarding this unknown metabolic crosstalk, our data pro- 
vide evidence supporting the hypothesis that raffinose repre- 
sents an important carbon source during fruit set and early 
development (Patrick and Offler, 1996; Nguyen-Quoc and 
Foyer, 2001), as its level was reduced in the source leaves of 
the transgenic plants. Raffinose is synthesized from UDP-Glc 
by subsequent reactions catalysed by UDP -glucose 4-epimer- 
ase (EC 5.1.3.2), galactinol synthase (EC 2.4.1.123), and raf- 
finose synthase (EC 2.4.1.82). The latter combines galactinol 
and sucrose to produce myoinositol and raffinose (Unda 
et al, 2012). In this regard, it has been demonstrated that, 
despite strong evidence for apoplastic unloading of sugars 
(Fridman et al, 2004; Hackel et al, 2006), a functional sym- 
plastic continuity of cell connections and numerous plas- 
modesmata between the phloem and storage parenchyma 
cells exists in young fruits. It is thus conceivable that carbon 
is partially unloaded via the symplast and, therefore, raffi- 
nose might represent an important carbon source. Virtually, 
all cell-wall precursors are derived from UDP-D-glucose 



(UDP-Glc) and GDP-D-mannose (GDP-Man) nucleo- 
tide sugars via a series of mechanistically related reactions 
that are cataysed by 4-epimerases, 3,5-epimerases, 4-reduc- 
tases, 4,6-dehydratases, dehydrogenases, and decarboxylases 
(Reiter, 2008). Thus, the reduction of GAUT4 activity could 
unbalance nucleotide sugar interconversion reactions lead- 
ing to a raffinose content reduction in transgenic lines once 
it is synthesized from UDP-Glc (Caffall and Mohnen, 2009). 
Moreover, the sucrose reduction trend observed in leaves of 
silenced plants re-enforces this hypothesis. 

Fruits of transgenic plants exhibited increased firmness. 
During tomato ripening, a decrease in the hydrostatic pres- 
sure (turgor) within fruit cells has been well documented to 
contribute to softening (Shackel et al, 1991). Presumably, this 
is due to a combination of events that lead to a dilution of 
the cell content and/or the direct loss of water from the fruit. 
These events may include: redistribution of solutes or solute 
leakage from the cell, breakdown of cell-wall structure, and 
transpiration. The elevation in apoplastic solute concentra- 
tions may occur as a result of increasing membrane leakiness 
and/or breakdown of complex cell-wall polysaccharides into 
smaller soluble components (Harker et al, 1997). Although 
membrane structure remains intact during fruit ripening, 
membrane conductivity (permeability) increases and thus 
there is a tendency for solutes to leak into the extracellular 
fluid, for the cell wall to hydrate, and for intercellular spaces 
to become water-soaked (Harker et al, 1997). However, the 
possible influence of solutes derived from cell-wall catabolism 
on the osmotic potential of the apoplast fluid has not been 
demonstrated. In this context, the higher firmness found in 
GylC/r^-silenced fruits could be explained by the following 
hypotheses. First, as a lower amount of ammonium oxalate- 
soluble pectin was observed, less intercellular solute accumu- 
lation might explain the higher cell turgor. Recently, Lunn 
et al. (2013) reported a similar phenotype where transgenic 
tomatoes with a blockage in pectin deposition showed a 
reduced amount of pectin together with an increase in fruit 
firmness. Secondly, the reduction in carbon influx observed 
in transgenic fruits could alter the increase in membrane 
conductivity, while the increment in free GalA contributes 
to water retention. Further experiments will be conducted to 
clarify the mechanisms underlying this phenomenon. 

Our results additionally bring insights to explain the impos- 
sibility of identifying gaut4 mutants in A. thaliana, described 
by Caffall et al. (2009). First, it has been reported that A. thal- 
iana seed mucilage is primarily composed of RG-I (Arsovski 
et al. , 2009); hence, a gaut4 mutation would compromise seed 
viability. Secondly, raffinose is an important source of carbon 
storage and desiccation tolerance in Arabidopsis seeds and 
consequently the alteration of nucleotide sugar interconver- 
tion reactions would compromise UDP-Glc production and 
ultimately raffinose biosynthesis (Zhou et al, 201 1). 

Interestingly, G/K/T^-silenced plants mimic the S. pennellii 
introgressed line 4-3-2 (IL4-3-2). This line exhibited a lower 
amount of GAUT4 mRNA (Fei et al, 2011), showed a fruit 
free GalA QTL, and exhibited dramatically reduced yield 
(Schauer et al. , 2006). Likewise, the transgenic plants analysed 
here displayed less pectin and an increase in free GalA in the 
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fruits, as well as a reduction in the harvest index. Therefore, 
the results shown in this work suggest that the GAUT4 gene 
is one of, or the sole, genetic determinant(s) of the S. pen- 
nellii introgressed line phenotype as we suggested previously 
(Bermudez et al. 2008). In the same genomic region of chro- 
mosome 4, Chai'b et al. (2006) mapped a QTL for fruit mea- 
liness and instrumental firmness. Furthermore, biochemical 
characterization revealed alterations in pectin composition as 
well as in cell size (Lahaye et al, 2012, 2013), in agreement 
with the results presented here. 

In conclusion, our results have presented characterization 
of the GAUTgsm family in tomato and showed that changes 
in the structure of cell-wall polysaccharides correlate with 
their spatial and temporal locations in developing plant tis- 
sues. Moreover, silencing of the GAUT4 gene altered pectin 
composition and affected plant growth and development. 
The observed phenotypes revealed physiological mechanisms, 
via cell-wall polysaccharides and raffinose metabolisms, that 
modulate plant resource allocation, resulting in an unprec- 
edented shift in source to sink carbon partioning in tomato. 

Supplementary data 

Supplementary data are available at JXB online. 

Supplementary Table SI. Primers used in qRT-PCR. 

Supplementary Table S2. Photosynthetic parameters meas- 
ured by gas exchange in tomato G^C/T^-silenced lines. 

Supplementary Table S3. Metabolite relative content in 
G^4t/7V-silenced lines. 

Supplementary Fig. SI. Heritability of the silenced 
phenotype. 

Supplementary Fig. S2. Leaf structural analyses of 
G^4t/7V-silenced plants. 
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